Electromagnons in the spin collinear state of a triangular lattice antiferromagnet 
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Terahertz time-domain spectroscopy was performed to directly probe the low-energy (1-5 meV) 
electrodynamics of triangular lattice antiferromagnets CuFei-^Ga^C^ (x = 0.00, 0.01, and 0.035). 
We discovered an electromagnon (electric-field-active magnon) excitation at 2.3 meV in the para- 
electric tt44 collinear magnetic phase, while this electromagnon vanishes in the ferroelectric he- 
limagnetic phase. Anti-correlation with noncollinear magnetism excludes the exchange-striction 
mechanism as the origin of dynamical magnetoelectric coupling, and hence evidences the observa- 
tion of spin-orbit coupling mediated electromagnon in the present compound. 

PACS numbers: 75.85.-K, 76.50.+g, 78.20.-e 



Magnetoelectric (ME) effect, i.e. electric (magnetic) 
control of magnetic (dielectric) properties, has long been 
an important issue in condensed matter physics [2} . While 
magnetic and dielectric properties usually show minimal 
coupling, recent discoveries of magnetically-induced fer- 
roelectricity in frustrated magnets have enabled unprece- 
dentedly large and versatile ME responses Jia et al. 
suggested that at least three microscopic ME coupling 
mechanisms can be considered: When a ligand ion is 
placed at the center of two magnetic ions, magnetically- 
induced local electric polarization Py is described as 



Pit = U(Si ■§ j ) + Aeij x {§i x Sj) 
+B[(eij ■ Si) Si - {cij ■ Sj)Sj], 



(1) 



where A and B are coupling coefficients, II is a vector 
unique to the underlying crystal structure, and eV,- an 
unit vector connecting two magnetic moments Si and Sj , 
respectively @. The first term represents the exchange 
striction, and is considered as the origin of ferroelectricity 
in some collinear antiferromagnets like Ca3CoMn06[7|. 
The second term comes from the inverse effect of 
Dzyaloshinskii-Moriya (DM) interaction^, which well 
explains the ferroelectricity in many noncollinear spiral 
magnets such as i?MnO 3 0, i] and Ni 3 V 2 8 [lf|. The 
third term arises from spin-dependent modulation of co- 
valency(hybridization) between metal d-state and ligand 
p-state[6|, while this term usually oscillates and cancels 
out within the crystal. Note that the second and third 
terms in Eq. (1) rely on the spin-orbit interaction, but 
the first term (exchange striction) does not[|j]. With any 
mechanism, a modification of magnetic structure, e.g. by 
external magnetic field H, leads to a significant change 
of induced electric polarization P. 

One important consequence of such a strong ME cou- 
pling is the appearance of a novel collective excita- 
tion called electromagnon (i.e. magnon driven by a.c. 
electric field E u ) in the dynamical regime. With de- 



tailed polarization analyses of absorption spectra, exis- 
tence of electromagnon (EM) excitation has been estab- 
lished for ferroelectric (FE) helimagnets -RMnC>3 [llj, ll2j , 
i?Mn 2 5 [l3j], and Ba 2 Mg2Fei 2 22 [lI|. Here > the most 
crucial is the microscopic origin of dynamical ME cou- 
pling, which is not necessarily identical to that of the 
magnetically-induced static P in the same compound. 
According to the inverse DM scheme, the EM excitation 
in FE helimagnets emerges as the rotational oscillation of 
spin-spiral plane and associated P- vector 15j . While this 
rotational mode should be active only with E u perpen- 
dicular to the spin-spiral plane, the observed selection 
rule for i?Mn 2 C>5 contradicted with this prediction 13, 
Il6j |. For i?MnC>3 and Ba 2 Mg 2 Fei 2 22 , the selection rules 
remain unchan ged even after the spin-flop transition un- 
der applied H 12|, [l4( . Latest theoretical studies sug- 
gested that the exchange striction mechanism can also 
host the EM activity in noncollinear magnets, but with 
the selection rule tied to the chemical lattice [17] . Impor- 
tantly, this exchange-striction-induced EM is inactive in 
the collinear magnetic phase, since the differential polar- 
ization 8 P^ oc Si ■ SSj becomes always zero (Si _L SSj). 
This latter model well reproduces the observed selection 
rules or absorption spectra of EM in all the three he- 
limagnetic compounds 14, 16- HI, whereas the firm ex- 
perimental evidence of spin-orbit coupling mediated EM 
excitation is still lacking. 

In this study, we report the experimental discovery 
of electromagnon excitation in the paraelectric collinear 
magnetic (ttii) phase of triangular lattice antiferro- 
magnets CuFei_^Gaa;0 2 . This EM mode was found 
to rather vanish in the FE helimagnetic phase. The 
anti-correlation between the electromagnon and the 
noncollinear magnetism excludes the exchange-striction 
mechanism as the origin of dynamical ME coupling, and 
hence suggests that the observed electromagnon is elec- 
trically activated by the spin-orbit coupling. 

CuFei_ x Ga 2: 2 crystalizes into a delafossite struc- 
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ture with centrosymmetric space group i?3m(Fig. 1(a)), 
which consists of the stacking of triangular lattices along 
the c-axis[l9j|. Each Fe 3+ (5 = 5/2) ion is surrounded 
by O 2- octahedra, and the magnetic frustration leads to 
several complex spin structures with propagation vector 
k = (?j 9) 3/2) [20(. For simplicity, we define the q- vector 
(q) as the in-plane component of magnetic propagation 
vector. CuFe02 is characterized by the collincar tt44- 
(CM4) magnetic ground state below 11K, with commen- 
surate q = 0.25 and spin direction along the c-axis(Fig. 
l(e))[21|. In Ga-doped specimen with x > 0.02, the 
CM4 phase is replaced by ferroelectric proper screw (NC) 
magnetic phase with incommensurate q ~ 0.202 [22h24| . 
where spin rotates within a plane perpendicular to the o- 
vector and P appears parallel to q | [110] (Fig. l(d))[25j. 
The ferroelectricity in the NC phase can be explained 
by neither exchange striction nor inverse D-M mecha- 
nism. Instead, Arima suggested that the third term in 
Eq. (1) can induce finite P \\ q component on the de- 
lafossite lattice with proper screw magnetic order [26}. 
This model predicts that the reversal of P-vector is cou- 
pled with the reversal of vector spin chirality (i.e. clock- 
wise or counter-clockwise manner of spin rotation), which 
was later confirmed by the polarized neutron scattering 
experiments (22|. The x - T phase diagram [2^ as well 
as the H - T phase diagram for the x = 0.01 specimen 
(by the present study) are summarized in Figs. 1(b) and 
(c). Note that the transition from CM4 into NC(FE) can 
also be induced by H applied along the c-axis[2(], H3- 
Starting from either of these magnetic ground states, in- 
crease of temperature (T) first induces partially disor- 
dered collinear (ICM1 or ICM2) magnetic phases (2lT [28 1 
and then produces a paramagnetic (PM) phase. All mag- 
netic phases other than NC(FE) are paraelectric[22l |23| . 

Single crystals of CuFei-zGa^O^ (x = 0.00, 0.01, and 
0.035) were grown by a floating zone method[19(. Mag- 
netization was measured with a SQUID magnetometer. 
Complex transmittance t is obtained by terahertz time- 
domain spectroscopy, and further converted into complex 
refractive index n — Jefi using the relationship 
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where e, //, d, u, and c represent the complex dielectric 
constant, complex magnetic permeability, sample thick- 
ness, frequency of light, and velocity of light, respec- 
tively. To numerically solve Eq. (2), we make the pre- 
cxponential factor approximate 4n/(n+ l) 2 by assuming 
H ~ 1 unless otherwise noted. The detail of THz time- 
domain spectroscopy is described in Appendix and Ref. 
12]. 

We first investigated the low-energy electrodynamics 
in the paraelectric CM4 collinear magnetic phase. Note 
that symmetry of triangular lattice allows the existence 
of three equivalent q || (110). The spectra observed in 
the present study reflect the contributions from all the 
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FIG. 1: (color online), (a) Crystal structure and (b) x - T 
magnetic phase diagram at H = for CuFei-^Ga^C^. (c) H 
- T magnetic phase diagram for the x — 0.01 specimen with 
static H applied parallel to the c-axis. Circles and squares 
are the data points obtained from the measurements of mag- 
netization with T- and //-increasing runs, respectively. Open 
(closed) symbols represent the data points determined in the 
present work (the previous work by Terada et al. 23]). (d) 
and (e) indicate magnetic structures of the NC(FE) and CM4 
phase, respectively. Dashed square in (e) represents the mag- 
netic unit cell in the CM4 phase. The directions of magnetic 
q-vector and electric polarization P are also indicated. 



three q-domains. Figures 2(a) and (b) indicate the real 
and imaginary part of eu spectra (Re[e/i] and Im[e/i]) 
with various polarization configurations for the x = 0.01 
specimen at 4.4 K, respectively. With E" || [110] and 
H w || [110], two resonance modes are observed at 1.2 and 
2.3 meV. Only the former one survives for E u \\ [001] 
and H u || [110], whereas only the latter one does for 
E u || [110] and H u \\ [001]. These results unveil that the 
excitation at 2.3 meV is an EM mode driven by E u \\ 
[110], while the one at 1.2 meV is a conventional magnon 
mode driven by H u || [110]. 

In the following, we focus on the behavior with the 
E u || [110] and H" \\ [110] configuration. To further ana- 
lyze the aforementioned e/i spectrum, the corresponding 
absorption coefficient a(= — 2(ln|t|)/d) and the decom- 
posed e and /x spectra are plotted in Figs. 3 (a)-(c), 
respectively. To discriminate the e- and /^-contributions 
to the e/i spectrum, we first assumed /t = 1 for hu > 2.0 
meV. The obtained e spectrum can be fitted well with 
the sum of two Lorentzian functions; a higher- frequency 
mode represents the lowest-lying optical phonon to give 
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FIG. 2: (color online). Real and imaginary parts of e/i spectra 
(Re [e/x] and Im [efj]) for the a; = 0.01 specimen measured at 
4.4 K with various light-polarization configurations. 



rise to the tail absorption observed below 5 meV. By 
substituting this e-fitting function into Eq. (2), the \x 
spectrum is deduced for Uui < 2.0 meV. It can be fitted 
as well with a single Lorentzian function. 

In Fig. 4(b), we show T-dependence of Im[e/i] spec- 
trum for the x = 0.01 specimen. Increase of temper- 
ature leads to broadening of two resonance peaks, and 
they become almost undiscernible in the ICM1 magnetic 
phase above 10 K. ICM1 and ICM2 are partially dis- 
ordered magnetic phases [U [28[, and may lose the spin 
correlation enough for magnons or EM excitations to be 
observed. The undoped x = 0.00 specimen also has 
two resonance modes at the same frequency in the CM4 
magnetic ground state (Fig. 4(a)), and shows similar 
T-dependence of Im[e/x] spectrum as observed for the 
x = 0.01 specimen. In contrast, the x = 0.035 specimen 
with the NC(FE) ground state shows no discernible peak 
structure in the whole temperature range (Fig. 4(c)). 
We also measured the Imfe^t] spectrum for the x = 0.035 
specimen with E || [001] and H || [110], but no peak 
structure was discerned. To summarize, the EM excita- 
tion driven by E u | [110] is active only in the paraelectric 
collinear CM4 magnetic phase, not in the ferroelectric NC 
helimagnetic phase. 

Next, we discuss the microscopic origin of these exci- 
tations. The spin-wave (SW) dispersion for CuFe02 has 
been investigated by a previous inelastic neutron diffrac- 
tion study [29(, and the analysis clarified the existence of 
two SW branches as reproduced in Fig. 3(g) [2(3]. In gen- 
eral, an antiferromagnetic resonance (AFMR) appears 
as the excitation of zone center mode at k = by H u 
perpendicular to the collinear spin direction(Fig. 3(e)). 
From this criterion, we concluded the excitation at 1.2 
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FIG. 3: (color online), (a)-(c) Absorption coefficient a, real 
and imaginary part of e and jj, spectra for the x = 0.01 spec- 
imen at 4.4 K with || [110] and H w || [llO]. Solid lines 
in (b) and (c) represent the fits with the sum of Lorentzian 
functions, (d) Spin structure of the CM4 magnetic ground 
state, (e) and (f) indicate the possible excitation modes cor- 
responding to the observed genuine magnon (AFMR) and elec- 
tromagnon(EM), respectively. In (f), Si and S3 rotate to the 
opposite direction of 52 and 54 within a plane perpendicular 
to q. (g) Spin-wave (magnon) dispersion of CuFe02 along the 
(h, h, 0) direction as suggested in Ref. [30| . 



meV driven by H u || [110] is AFMR on the lower SW 
branch. In contrast, the excitation energy of the ob- 
served EM (~ 2.3 meV) agrees with that of the zone- 
center mode on the upper SW branch. 

So far, the most successful scheme to explain the dy- 
namical ME coupling is the exchange striction. How- 
ever, this mechanism is inactive in the collinear spin 
system like the present CM4 phase, since the relation- 
ship SPij oc Si ■ SSj = always holds (l7L liH ]. This 



strongly suggests the relevance of spin-orbit coupling to 
the present EM mode; the relatively weak peak intensity 
in the Im[e] spectrum (one order of magnitude smaller 
than that of DyMn03[12|) also supports this scenario. 
Since the static P in NC(FE) is induced by the proper 
screw magnetic order through the spin-orbit-interaction 
mediated modulation of Fe 3d - O 2p hybridization 0, [2^], 
we may anticipate the analogous origin for the presently 
observed dynamical ME coupling in CM4. For example, 
the magnetic excitation as depicted in Fig. 3(f) can dy- 
namically generate a proper-screw-like spin texture with 
a finite spin chirality, which is expected to induce non- 
zero electric dipole along the SP || q || [110] direction. 
This mode should be active only with E u || [110], which 
is consistent with the experimental results. The disap- 
pearance of EM mode in the NC(FE) phase may reflect 
the alteration of magnetic symmetry or Brillouin-zone 
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FIG. 4: (color online). T-dependence of Im[e^t] spectra mea- 
sured with E u || [110] and H" || [110] for (a) x = 0.00, (b) 
x = 0.01, and (c) x = 0.035 specimens, respectively, (d) 
Im[e] spectra for the x — 0.01 specimen measured at 6 K with 
W || [110] and H" || [001] in various magnitudes of static H 
(i/ DC ) applied along the [001] direction, //-dependence of ob- 
served EM peak positions (circle), as well as the development 
of resonance modes previously reported by the ESR study for 
the x — 0.00 specimen [§2] (square), are plotted in (e). 



folding, but the detail is left to be clarified. 

We further investigated the development of the EM 
mode in static H applied along the [001] direction. Figure 
4(d) indicates the //-dependence of Im[e] spectrum mea- 
sured at 6 K for the x = 0.01 specimen with E u || [110] 
and H u | [001], where only the EM excitation can be ob- 
served. As H increases, the EM mode is found to split lin- 
early with H and form two peak structures. This reflects 
the //-linear splitting of SW branches, which is generally 
expected for collinear antiferromagnets with H applied 
parallel to the magnetic easy axis. Since the spectral 
shape of the two-magnon excitation should be indepen- 
dent of -ff[3l|, this ensures that the present EM mode is 
excited by the one-magnon process. The observed evolu- 
tion of EM peak positions under applied H is summarized 
in Fig. 4(e). Note that similar ii-dependence of reso- 
nance modes has been reported by Fukuda et al. from 
the ESR study for the x = 0.00 specimen [32|, while they 
conventionally assigned these modes to AFMR driven 
by H u . Our present results imply that the resonance 
modes found in the previous ESR study [3^] is primar- 
ily driven by i?"-componcnt of incident microwave. The 
peak structure observed in Im[e] spectrum becomes al- 
most invisible after the transition from CM4 into NC(FE) 
at 6.3 T, which confirms the inactivity of EM mode in 



the latter NC phase. 

In conclusion, we have experimentally revealed 
the clcctromagnon excitation in the paraclcctric 
collinear magnetic phase of triangular lattice antifcrro- 
magnet CuFei_ a; Ga a ;02. This mode was found to vanish 
in the ferroelectric helimagnetic phase. These facts prove 
that neither ferroelectricity nor noncollinear magnetism 
is a necessary condition for the appearance of electro- 
magnon excitation, while the existing theories on electro- 
magnon have focused on noncollinear magnets like heli- 
magnets. The anti-correlation between the noncollinear 
magnetism and the emergence of electromagnon excludes 
the exchange striction mechanism as the origin of dynam- 
ical ME coupling: The electric activity of the magnon in 
this compound is ascribed to the modulation of the p- 
d hybridization at the spin-twisted excited state via the 
spin-orbit interaction. Our discovery suggests that sim- 
ilar electromagnon modes will be observable in a wide 
range of paraelectric collinear magnets. 



APPENDIX 

In this appendix, we provide the detail of the tera- 
hertz time-domain spectroscopy and discuss the validity 
of our transformation procedure among different optical 
constants. 

In the present study, terahertz time-domain spec- 
troscopy (THz-TDS) was performed in the transmission 
geometr y. T he detailed experimental setup is described 
in Ref. [l2j. We employed the photo-conducting anten- 
nas made of low-temperature-grown GaAs as emitters 
and detectors, to access the energy down to 1 meV. The 
direction of the light-polarization was carefully set paral- 
lel to the crystallographic axis using a wire grid polarizer. 

In THz-TDS, wave form of irradiated pulse electric 
field (ranging within a few picoseconds) is directly mea- 
sured in time domain with and without the specimen. 
They are converted into the frequency domain via the 
fast Fourier transformation (FFT), and the spectra of 
complex transmission constant (t) is deduced using the 
relationship t = E" ample / E" e{ . Here, we can avoid the ef- 
fect of multiple reflections and resultant interference by 
restricting the time range of FFT. Obtained t is further 
transformed into complex refractive index n — ^JeJL us- 
ing Eq. (2). Since Eq. (2) is not explicitly solvable, we 
made its pre-exponential factor approximate 4n/(n + l) 2 
by assuming [i ~ 1 unless otherwise noted. 

This approximation hardly affects the obtained n (or 
e/i) spectrum, as exemplified in the following by the spec- 
tra for the x = 0.01 specimen with E u || [110] and 
H 1 ^ || [110] at 4.4 K. A solid line in Fig. 5 indicates the 
imaginary part of e/i-spectrum deduced with the afore- 
mentioned (i ~ 1 approximation. Its decomposed e- and 
/i-spectra, which are derived with the method described 
in the main text, are also shown in Figs. 3 (b) and (c). 
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FIG. 5: (color online). Imaginary part of e/j, spectra for the 
x = 0.01 specimen measured at 4.4 K with E" || [110] and 
H u || [110] in zero magnetic field. For a solid (dashed) line, 
fi = 1.0 (/x = 1.05 + 0.08z) is assumed for the pre-exponential 
factor in Eq. (2), respectively. Note that the solid line repre- 
sents the same data as shown in Fig. 2(b). 



We find that the relationship \fx — 1| < 0.1 always holds, 
which justifies the presently assumed /i ~ 1 condition. 
To further check the validity of this approximation, we 
deduced e/i spectrum assuming /i = 1.05 + 0.08i (the 
maximum (i- value taken from Figs. 3 (b) and (c)) for 
the pre-exponential factor in Eq. (2). The result is plot- 
ted in Fig. 5 as a dashed line, which is consistent with 
the one calculated with the original /i = 1 assumption 
(the solid line). 
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